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Abstract. We present a novel model to simulate real social networks of complex interactions, based in a granular 
system of colliding particles (agents). The network is build by keeping track of the collisions and evolves in time with 
correlations which emerge due to the mobility of the agents. Therefore, statistical features are a consequence only of 
local collisions among its individual agents. Agent dynamics is realized by an event-driven algorithm of collisions 
where energy is gained as opposed to granular systems which have dissipation. The model reproduces empirical data 
from networks of sexual interactions, not previously obtained with other approaches. 
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PACS. 89.65.-s Social and economic systems - 89.75 .Fb Structures and organization in complex systems - 
Networks and genealogical trees - 89. 75. Da Systems obeying scaling laws 
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1 Introduction 

A social network is a set of people, each of whom is acquainted 
with some subset of the others. In such a network the nodes (or 
vertices) represent people joined by edges denoting acquain- 
tance or collaboration. Empirical data of social networks in- 
clude networks of scientific collaboration! l|i of film actor co- 
llaborations^!, friendship networks! 3 1 among some others! 4 1. 
One kind of social network is the network of sexual contacts! 5 
161718 1. where connections link those persons (agents) that have 
had sexual contact with each other. The empirical investigation 
of such networks are of great interest because, e.g. the topo- 
logical features of sexual partners distributions help to explain 
why persons can have the same number of sexual partners and 
yet being at distinct risk levels of contracting HIV 1 9 1 . 

The simplest way to characterize the influence of each in- 
dividual on the network is through its degree k, the number 
of other persons to whom the individual is connected. Sexual 
contact networks are usually addressed as an example of scale- 
free networks I6I7I81 1 01 . because its properly of having a tail 
in its degree distribution, which is well fitted by a power-law 
P(k) ~ k^ 1 ', with an exponent 7 between 2 and 3. However, 
another characteristic feature, not taken into account, is that 
the small fc-region, comprehending the small-k values varies 
slowly with k, deviating from the power-law. Moreover, the 
size of the small-A; region also increases in time, yielding rather 
different distributions when considering the number of partners 
during a one year period or during the entire life, e.g. for entire- 
life sexual contacts, the degree distribution shows that at least 
half of the nodes have degree in the small-/c region 15 161 . A 
model predicting all these different distributions shapes for dif- 
ferent time spans is of crucial interest, because the transmission 
of diseases occur during the growth mechanism of the network. 
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Fig. 1. Snapshots of the growing network of collisions in a low-density 
gas with N = 100, for (a) n = 0.02iV, (b) n = 0.15/Y, (c) n = N. 
In the on-line version two colors represent the two sexual genders and 
larger symbols emphasize those which belong to the network (linked 
agents). 



One of the main difficulties for validating a model of sex- 
ual interactions is that typical network studies of sexual con- 
tacts involve the circulation of surveys, or anonymous ques- 
tionnaires, and only the number of sexual partners of each in- 
terviewed person is known, not being possible to obtain infor- 
mation about the entire network, in order to calculate degree 
correlations, closed paths (cycles), or average distance between 
nodes. 

In this work we propose a model of mobile agents in two 
dimensions from which the network is build by keeping track 
of the collisions between agents, representing the interactions 
among them. In this way, the connections are a result not of 
some a priori knowledge about the network structure but of 
some local dynamics of the agents from which the complex 
networks emerge. Below, we show that this model is suitable 
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to reproduce sexual contact networks with degree distribution 
evolving in time, and we validate the model using contact trac- 
ing studies from health laboratories, where the entire contact 
network is known. In this way, we are able to compare the num- 
ber of cycles and average shortest path between nodes as well 
as compare the results with the ones obtained with Barabasi- 
Albert scale-free networks 1 1 1 1, which are well-known models, 
accepted for sexual networks I6I7I8I . We start in Sec.|2]by de- 
scribing the model of mobile agents, and in Sec. [5] we apply 
it to reproduce the statistical features of empirical networks of 
sexual contacts. Discussion and conclusions are given in Sec. [4] 

2 Model of mobile agents for sexual 
interactions 

The model introduced below is a sort of a granular system| 12), 
where N particles with small diameter d represent agents ran- 
domly distributed in a two-dimensional system of linear size 
L ^/N d (low density) and the basic ingredients are an in- 
crease of velocity when collisions produce sexual contacts, two 
genders for the agents (male and female), and n/N, the frac- 
tion of agents that belong to the network, which constitutes an 
implicit parameter for the resulting topology of the evolving 
the network. 

The system has periodic boundary conditions and is ini- 
tialized as follows: all agents have a randomly chosen gender, 
position and moving direction with the same velocity modulus 
| wo I- We mark one agent from which the network will be con- 
structed. When the marked agent collides for the first time with 
another one of the opposite gender, the corresponding collision 
is taken as the first connection of our network and its colliding 
partner is marked as the second agent of the network (Fig.[^). 
Through time, more and more collisions occur, increasing the 
size n of the network (Fig. [Q) and[0;) till eventually all the 
agents composing the system are connected. 

Collisions between two agents take place whenever then- 
distance is equal to their diameter and the collision process 
is based on an event-driven algorithm, i.e. the simulation pro- 
gresses by means of a time ordered sequence of collision events 
and between collisions each agent follows a ballistic trajecto- 
ry! 13 1. Since sexual interactions rely on the sociological ob- 
servation! 14 1 that individuals with a larger number of partners 
are more likely to get new partners, we choose a collision rule 
where the velocity of each agent increases with the number k of 
sexual partners. The larger the velocity one agent has the more 
likely it is to collide. Moreover, contrary to collision interac- 
tions where velocity direction is completely deterministic! 15 1, 
here the moving directions after collisions are randomly se- 
lected, since in general, sexual interactions do not determine 
the direction towards which each agent will be moving after- 
wards. Therefore, momentum is not conserved. 

Regarding these observations our collision rule for sexual 
interactions reads 

v{h) = {k? + \v (i)\)u, (1) 

where fc, is the total number of sexual partners of agent i, expo- 
nent a is a real positive parameter, u) — (e x cos 6 + e y sin 8) 




Fig. 2. (a) Cumulative distribution P C um (k) of the number k of part- 
ners among agents, when considering type-(i) and -(ii) interactions 
(see text) for n = 10 3 (circles), n = 10 4 (squares), n = 5 x 10 4 
(diamonds) and n = 10° (triangles). For the same parameter values 
(b) shows a pure scale-free distribution, obtained when only type-(ii) 
interactions form links. The solid line indicates the slope 7 = 3 of the 
scale-free distribution. Here a = 1 and N = 320 x 320. 

with 8 a random angle and e x and e y are unit vectors. Col- 
lisions which do not correspond to sexual interactions only 
change the direction of motion. 

Collisions corresponding to sexual interactions, i.e. with a 
velocity update as in Eq. 0, are the only ones which produce 
links, and occur in two possible situations: (i) between two 
agents which already belong to the network, i.e. between two 
sexually initiated agents and (ii) when one of such agents finds 
a non-connected (sexually non-initiated) agent. For simplicity, 
we do not take into account sexual interactions between two 
non-connected agents, and therefore our network is connected 
(see the discussion in Sec.|4). 

When interactions of type (i) and (ii) occur, both the dis- 
tribution tail and the small-fc region are observed, as shown by 
the cumulative distribution P cum {k) in Fig.|2^. Here, we use a 
system of N = 320 x 320 agents with p = 0.02, a = 1 and dis- 
tributions are plotted for different stages of the network growth, 
namely n = 10 3 , n = 10 4 , n = 5 x 10 4 and n = 10 5 - N. 
As one sees, the exponent of the power-law tail and the tran- 
sition between the tail and the small-fc region increase during 
the growth process. These features appear due to the fact that at 
later stages most of the collisions occur between already con- 
nected agents. Consequently, the average number of partners 
increases as well. 

If one considers only type-(ii) sexual contacts, the system 
reproduces a stationary scale-free network, as shown in Fig.|2j3. 
In this case the average number of partners, defined as| 16 1 
(fc) = fc TOJ :„(7 — l)/(7 — 2) with k m i n the minimum num- 
ber of partners, is always 2 (k m i n — 1 and 7 = 3). As we show 
below, while empirical data of sexual contacts over large peri- 
ods have distributions like the ones for regime (i)+(ii), data for 
shorter periods (1 — 10 years) are scale-free (only (ii)). 
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Fig. 3. Cumulative distributions, when varying (a) a parameter s of 
selectivity which interpolates between Figs|2^ and|2p (see text) for 
q = 1 and (b) the exponent a in the update velocity rule, Eq. Q, 
for s — 0. Here the same stage of growth is considered, namely n = 
5 x 10 4 = 0.5/V. 



With our model one can easily interpolate between both 
interaction regimes, (i)+(ii) and (ii), by introducing a parameter 
s of 'selectivity', defined as the probability that sexual initiated 
agents in case of collision with another initiated agent, have 
no sexual contact. Physically, this selectivity accounts for the 
intrinsic ability that a node has to select from all its contacts 
(collisions) the ones which are sexual. These intrinsic abilities 
were already used in other contexts, e.g. as a new mechanism 
leading to scale-free networks in cases where the power-law 
degree distribution is neither related to dynamical properties 
nor to preferential attachmentl 17 1. For s = one obtains the 
two regions illustrated in Fig. [2^, namely the small- k region 
and the power-law tail, while for s = 1 one obtains the pure 
scale-free topology illustrated in Fig. |2j). In Fig. fJt, we show 
the crossover between these two regimes. 

The shape of the cumulative distributions is also sensible to 
the exponent a in the update velocity rule, Eq. {0, as shown in 
Fig-Eb- While for small values of a < 1 one gets an exponential- 
like distribution, for a > 1.4 the distribution shows that a 
few nodes make most of the connections. Henceforth, we fix 
a = 1.2. 

Having described the model of mobile agents we proceed 
in the next section of a specific application, i.e. modeling em- 
pirical networks of sexual contacts. 



3 Reproducing networks of sexual contacts 

In this Section we will show that, by properly choosing the pa- 
rameter values in our model, one can reproduce real data dis- 
tributions of sexual contact networks. In Fig.|4]the cumulative 
distributions of a real contact network|6 1 are shown for females 
(circles) and males (triangles) separately, based on empirical 
data from 2810 persons in a Swedish survey of sexual behav- 
ior. The solid lines in Fig.|4]indicate the simulated distributions. 
The simulated power-law tails have exponents 7 ? „ = 2.4 and 
7/ = 4.0 for males and females respectively, compared with 
the empirical data -f m = 2.6 ± 0.3 and 7/ = 3.1 ± 0.3[6|. 



To stress that, while the power-law tails are also well fitted by 
distributions obtained with scale-free networks (dashed lines 
in Fig.|4j, these distributions have a minimum number of con- 
nections (partners) of fc m m = 5 for females and k m ; m = 7 
for males, contrary to the real value k m ; m = 1 also repro- 
duced with our agent model. In fact, the model of mobile agents 
takes into account not only the power-law tail of these distribu- 
tions, but also the small-fc region which comprehends the sig- 
nificant amount of individuals having only a few sexual part- 
ners (k > 1). 

Is important to note that in order to reproduce the differ- 
ence in the exponents is necessary to have 58% females and 
42% males, which is far from the expected difference among 
number of females and males in typical human populations, 
with ratios of females:males of the order of 1.1. The difference 
in the exponents of the distributions tails for males and females 
separately, present in the data of sexual surveys, has generated 
much controversy and is often considered due to a bias either 
of sampling or honest reporting (see Ref. |5| and references 
therein). The exponents j m and 7/ on a bipartite network are 
expected to be nearly the same when the percentage of males 
and females are similar, as shown in Fig. [5] In each plot five 
different realizations are shown for males (solid lines) and fe- 
males (dashed lines). In Fig.|5t, we plot the results for the same 
conditions as in Fig. E|(58% females and 42% males). Taking 
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Fig. 4. Cumulative distribution of sexual partners in a network of het- 
erosexual contacts extracted from Ref. [6], where male (triangles) and 
females (circles) distributions are plotted separately, with a total of 
2810 persons. Solid lines indicate the simulations when plotting the 
distributions at the same stage n — 0.27V, starting with a population 
composed by 58% of females and 42% of males. Here N = 10 , 
s = 
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Fig. 5. Distributions of sexual partners in a network of heterosexual 
contacts using different amounts of females: (a) 58% as in Fig.|4] (b) 
52%, (c) 64% and (d) 70%. In each plot solid and dashed lines in- 
dicate the cumulative distributions of males and females respectively, 
for five different realizations. Clearly, the exponent of the power-law 
tail of the distributions decreases when the percentage of females or 
males increases (see dotted lines). Same conditions as in Fig.[4]were 
used. 



, f . *- -# Scale Free 
<#t> o o Agents 



_ ^ a- - pt O^X 

jct ' 

-o 



n 



io 4 



10' 



Fig. 6. Comparing the average number of partners for scale-free net- 
works (stars) and the agent model (circles). For the scale-free network 
kmin = 4 and for agent model N = 520 x 520, s — 0. 





Fig. 7. Sketch of two real sexual contact networks having (a) only 
heterosexual contacts (TV = 82 nodes and L — 84 connections) and 
(b) homosexual contacts (TV = 250 nodes and L — 266 connections). 
While in the homosexual network triangles and squares appear, in the 
heterosexual network triangles are absent (see Table0. 
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Fig. 8. (a) Cumulative degree distribution of a homosexual contact 
network| 18 1 with n = 250 (triangles), (b) Cumulative degree distri- 
bution of a heterosexual contact network| 19 1 with n = 82 (triangles). 
Each case is compared with the average degree distribution over 20 
iterations, for the B A scale-free model (dashed line) with k m in = 1, 
and k m in — 2 and with our agent model (solid line) with s — 0.7. (c) 
Cluster coefficient for the homosexual network empirical data (trian- 
gles), the agent model (solid line) and the sf model k m i n = 2. The 
scale-free k m i n — 1 yields C(k) = (not shown). 



the average of the five curves for each gender yields the curves 
shown in Fig. [4] Figures |5j3-d show the distributions obtained 
for other percentages illustrating that when decreasing the dif- 
ference in the ratio females:males difference in the exponents 
dissappear. A characteristic feature of our model is that the av- 
erage number (k) of partners increases as the network grows, 
which is natural characteristic expected to occur in real sexual 
networks according to the observed differences in the shape 
of the degree distribution for yearly and entire-life reports of 
number of sexual partners |6 7|. This feature is not observed 
in scale-free networks, as illustrated in Fig. [6] Of course, that 
this growth also indicates non-stationary regimes, where (k) di- 
verges with the network growth. In the next section we explain 
how to overcome this shortcoming. 



We compare the model, with two empirical networks of 
sexual contacts. One network is obtained from an empirical 
data set, composed solely by heterosexual contacts among n = 
82 nodes, extracted at the Cadham Provincial Laboratory (Man- 
itoba, Canada) and is a 6-month block data 1 1 9 1 between Novem- 
ber 1997 and May 1998 (Figure^ sketches this network). The 
other data set is the largest cluster with n — 250 nodes in the 
records of a contact tracing study 1 18 1, from 1985 to 1999, for 
HIV tests in Colorado Springs (USA), where most of the reg- 
istered contacts were homosexual (see Figure^)). 

Figures EJa)-(b) show the cumulative distribution of the 
number of sexual partners for each of the empirical networks. 
For both cases the agent model and scale-free networks with 
kmin = 1 can reproduce the distribution of the number of 
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83.63 





1.45 





Homosexual 
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250 


287.03 


8.23 


10.52 


0.02302 


Heterosexual 
(Scale-free) 


82 


162 





159.72 





Homosexual 
(Scale-free) 


250 


498 


45.28 


256.79 


0.08170 



Table 1. Clustering coefficients and cycles in two real networks of 
sexual contacts (top), one where all contacts are heterosexual and an- 
other with homosexual contacts. In each case one indicates the values 
of the number N of nodes, the number L of connections, the number 
T of triangles, the number Q of squares and the average clustering 
coefficient C. The values of these quantities are also indicated for net- 
works constructed with the agent model and for scale-free networks 
with kmin ~ 2, note that for k m in = 1, L=81 and L=249, respec- 
tively and there are not cycles. 



partners. However, the agent model with s = 0.7 reproduces, 
as well, the clustering coefficient distribution that we measure 
from the empirical network. 

The clustering coefficient C(i) of one agent is defined|2| 
as the total number of triangular loops of connections passing 
through it divided by the total number of connections fcj . Aver- 
aging C(i) over all nodes with k{ neighbors yields the cluster- 
ing coefficient distribution C(k). While for the scale-free graph 
which better reproduces these empirical data, the clustering co- 
efficient is zero, our agent model yields a distribution which 
resembles the one observed in the real network (Fig.|SJ;). This 
feature is due to the co-existence of a tree-like substructure and 
closed paths (see Figs. 03). 

For both heterosexual and homosexual networks of sexual 
contacts, the model of mobile agents reproduces other impor- 
tant statistical features, namely the average clustering coeffi- 
cient C and the number of loops of a given order. Table [2 in- 
dicates the number T of triangles (loops composed by three 
edges), the number Q of squares (loops with four edges) and 
the average clustering coefficients C given by|2| the average 
of C(i) over the entire network. 

When using the agent model with the same number N of 
nodes as in the real networks we obtain similar results for L, T, 
Q and C, as shown in Table^middle), where values represent 
averages over samples of 100 realizations. For the heterosex- 
ual network there are no triangles due to the bipartite nature of 
the network. At the bottom of Table ^ we also show the val- 
ues obtained with scale-free networks, for both cases of one 
and two genders, whose minimum number of connections was 
chosen to be k m i n = 2, for which the clustering coefficient dis- 
tributions are as close as possible from the distributions of the 
real networks. Clearly, the agent model not only yields cluster- 
ing coefficient values much closer to the ones measured in the 
empirical data, but also does not show the formation of huge 



amounts of loops (triangles and squares), a feature of scale-free 
networks which is not observed in empirical data. 



4 Discussion and conclusions 

In this paper we presented a new model for networks of com- 
plex interactions, based on a granular system of mobile agents 
whose collision dynamics is governed by an efficient event- 
driven algorithm and generate the links (contacts) between a- 
gents. As a specific application, we showed that the dynamical 
rules for interactions in sexual networks can be written as a 
velocity update rule which is a function of a power a of the 
previous contacts of each colliding agent. For suitable values 
of a and selectivity s, the model not only reproduces empirical 
data of networks of sexual contacts but also generates networks 
with similar topological features as the real ones, a fact that is 
not observed when using standard scale-free networks of static 
nodes. 

Furthermore, our model predicts that the growth mecha- 
nism of sexual networks is not purely scale-free, due to inter- 
actions among internal agents, having a mean number of part- 
ners which increases in time. This should influence the predic- 
tions from models of spreading of infections 1 16 1. Our agent 
model offers a realistic approach to study the emergence of 
complex networks of interactions in real systems, using only 
local information for each agent, and may be well suited to 
study networks in sociophysics, biophysics and chemical re- 
actions, where interactions depend on specific local dynamical 
behavior of the elementary agents composing the network. 

While given promising results the model may be improved 
in two particular aspects. First, it should enable the conver- 
gence towards a stationary regime with a growth process start- 
ing with all possible collisions instead of one particular agents 
from which the network is constructed. Second, the depen- 
dence of the above results on the velocity rule in Eq. Q should 
be studied in detail, namely for the case of constant velocity 
{a = 0). Preliminary results have shown that the stationary 
regime is easily obtained with the model above by introducing 
a simple aging scheme, while by varying the parameter a one is 
able to reproduce other non-trivial degree distributions. More- 
over, we introduced the selectivity parameter s to select from 
all possibles social interactions (collisions) the ones which are 
of sexual nature. Without introducing this selectivity, the model 
of mobile agents is able to reproduce other social networks 
of acquaintances. These and other questions will be addressed 
elsewhere|20|. 
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